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1 Introduction 



. In a recent paper we proposed a realistic interpretation of the Schrodinger 



' equation for density matrices, in which the difference between the position argu- 

ments of the density matrix is considered as an objective extra space dimension. 
I> '. In the case of a free particle, where the potential V (x) vanishes, we found so- 

I lutions which are perfectly localized both in position space and in momentum 

' space; these solutions behave exactly as non-relativistic point-like particles mov- 

. ing at constant speed, with the correct values for all observable quantities. In 

Oh! the general case, where V (x) ^ 0, we pointed out that the natural frequencies 

of the stationary states in the density matrix representation correspond to the 
■ difference between two energy levels in the original quantum system; since the 

^ , "jumps" between energy levels are observable, while the individual energy levels 

are not, we trivially deduced that the observable natural frequencies are the same 
in both representations. 

^ , In the first paper the non-relativistic case was treated mainly as an introduc- 

I tion to the relativistic case; in this second paper we will study it in more detail, 

examining the correspondence between our new representation and the standard 
representation of non-relativistic quantum mechanics, based on the Schrodinger 
equation for pure states. In Section ^ we will first summarize the main results of 
the previous paper, and then we will give a qualitative definition of a quantum 
trajectory (describing an individual physical system), as opposed to a statistical 
mixture (describing an ensemble of physical systems); as a consequence, we will 
reject two basic principles of standard quantum mechanics, namely the superposi- 
tion principle and the wave function collapse. In Section |3| we will show the results 
obtained by applying our approach to simple potentials: by examining the linear 
harmonic oscillator we will shed a new light on the supposed equivalence between 
energy eigenvectors and stationary states; besides, the delta barrier potential will 
provide a semi-classical explanation of the tunnel effect, in which an essential fea- 
ture is the introduction of internal degrees of freedom even for spinless particles. 
Finally in Section ^ we will present our conclusions. 
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2 The density matrix representation and the superpo- 
sition principle 

The standard quantum description of a non-relativistic particle moving in one 
space dimension is based on a wave function ^ (x), whose time evolution is given 
by the Schrodinger equation: 

We will refer to equation as to the Schrodinger equation for pure states. The 
time evolution generated by @) belongs to the group of unitary transformations; 
a new representation of this group, the density matrix representation, is obtained 
by means of the fundamental relation 

V{x,y)=i;{x)r{y) (2) 

If we apply an arbitrary unitary transformation to the pure state ip{x), the re- 
lation (|2|) enables us to obtain the same transformation in the density matrix 
representation. Specifically, if we consider the hamiltonian operator H, the mo- 
mentum operator P and the position operator Q as generators respectively of time 
translations, space translations and momentum translations, we easily obtain the 
following expressions for the same generators in the new representation: 

Qip = {x -y)ip (3) 



,dx dy , 

The time evolution in the density matrix representation is then simply 

In the usual interpretation, equation (^) describes the time evolution of mixed 
states, and therefore its solutions have only a statistical meaning: they are simply 
useful tools for computing probability amplitudes and expectation values. On the 
contrary, in our approach the solutions of equation ^ are considered as real ob- 
jective fields describing individual physical systems. To avoid misunderstandings, 
in the rest of this paper we will use the expression "density matrix" when we mean 
the usual statistical interpretation and we will use the new expression "quantum 
matrix" when we mean a real objective field. 

A quantum matrix ip (x, y) depends upon two position coordinates x and y; if 
we want to consider (p{x,y) as a real objective field then we must introduce an 
objective extra space dimension. Therefore, we will define a new pair of position 
coordinates ^ 

xs = ^{x + y) XD = x-y (5) 
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and we will interpret xs as the "physical" position coordinate, while xd will 
be an "auxiliary" position coordinate, having observable effects only around the 
point X£) = 0; both xs and xd, however, will be considered as objective position 
coordinates (note that in the definition of xd we have inverted the sign convention 
with respect to our first paper [jl8|). In the rest of this paper we will write loosely 
if {x, y) or if {xs,xd), meaning the same field expressed in two different coordinate 
systems. 

The transition from position space to momentum space is obtained by means 
of the two-dimensional Fourier transform: 

'^ikx,ky) = ^ f [ if (x, y) e-^^^+'^y^dx dy (6) 



Again we define a new pair of momentum coordinates: 

ks = ^ (kx + ky) kD = kx- ky (7) 

and again ks will be the "physical" momentum coordinate, while ko will be an 
"auxiliary" momentum coordinate, having observable effects only around the point 

kD = 0. 

Since we have two position coordinates and two momentum coordinates, we 
naturally can define two position operators and two momentum operators. The 
operators already defined in (^) may be rewritten as: 

QDip = XD(p FD(p = -ih- — (8) 

dxs 

In addition, we introduce the new operators and by means of the following 
definition: 

d 

OXD 

In momentum space the operators Pd and P5 are simply given by: 

Pj^<^ = hkD'^ Fs'^ = ^ks<^> (10) 

It is very important here to note that Q5 and P5 are commuting operators, i.e. 
[Q5, P5] = 0; this property is fundamental for our approach, because it means that 
they have common eigenvectors. These eigenvectors have the form 



or, in momentum space. 



<^ = 5(x5-xo)e^'=°"^ (11) 



^ = 5 {ks - ko) e-''''^'"' (12) 



Since xs and ks are the "physical" position and momentum coordinates, it is clear 
that the field defined in (|ll|) and ([l^ ) is perfectly localized in position space at 
Xs = xq and in momentum space at ps = Po = ^^o- Therefore in the density 
matrix representation the Heisenberg uncertainty principle is not true and we 
conclude, as in our previous paper, that the Heisenberg principle has no ontological 
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meaning; rather, it is just a shortcoming of the standard quantum formahsm based 



on pure states. A similar point of view is expressed by Olavo |16|. 

While the two operators Q^) and are interesting just because they generate 
translations respectively in momentum space and in position space, the two opera- 
tors Q5 and P5 are more strictly related to the measurement of the corresponding 
physical quantities, as we will now see. To define the observable quantities in the 
density matrix representation, we start from the expressions for the mean values 
in the original representation and then exploit the fundamental relation (^); we 
obtain: 



Q 



Z£l=0 



(13) 



P 



(Ps 95) „dx5 



(14) 



E 



2m 



p| + y(Q, 



(15) 



a: 15=0 



The above definitions confirm our interpretation of as an "auxiliary" position 
coordinate, having observable effects only around xd = 0. The same quantities 
may be written in momentum space, obtaining analogous expressions; for instance, 
the momentum P may be written as 



(P5$) ^ dks 

ko=0 



(16) 



Since we do not assign a statistical meaning to the quantum matrix ip, we will 
not interpret the expressions (^3|)-(^) as mean values. On the contrary, we will 
consider (|l^) as the center of mass of the system, while (p^) and ( [l5| ) will be 
respectively the total momentum and the total energy. 

If we Fourier transform the quantum matrix if with respect to alone, we 
obtain the well known Moyal-Wigner transformation: 



F {xs,ps) 



(f {xs,xd) e 



dxD 



(17) 



first introduced by Wigner ||2^ in 1932. The Wigner function F{xs,ps) has 
some valuable properties: the common eigenvectors of the operators Q5 and 
Fs are represented simply by products of delta functions, i.e. F{xs,ps) = 
6 {xs — xq) 5 {ps — Po); besides, those operators which depend only on Qg and 
P5 become c-numbers in the Wigner representation; finally, if G is an ob- 
servable quantity depending on position and momentum through the function 
g{x,p) = gi{x) + g2(,p), then in the Wigner representation we can write the 
expression 

G = J J g{xs,Ps)F{xs,Ps)dxsdps (18) 

i.e. F {xs,Ps) weighs the function g{xs,Ps) over the whole phase space {xs,Ps), 
the application of ( p^ ) to the quantities defined in (|l3|)-(^) is straightforward. 
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Thus the Wigner function is a useful tool for establishing relations between the 
quantum description and the classical description based on phase space; indeed, 
it provides the basis for Moyal's deformation quantization |]l^, which is an au- 
tonomous formulation of quantum mechanics, alternative to the more familiar 
Hilbert space and path integral quantizations. 

However, in the usual statistical interpretation of the density matrix, expres- 



sion (18) yields the mean value of the quantity G, and therefore the Wigner 
function seems to play the role of a joint probability distribution in phase space. 
But even if the function (p (x, y) satisfies all the properties required for being a 
well defined density matrix, it is in any case possible for the Wigner function to be 
negative in some phase space region. This is the well known "negativity problem" , 
which prevents a complete analogy between the classical description in phase space 
and the quantum description in the Wigner representation. On the contrary, if 
we consider the quantum matrix ip (x, y) as a real objective field, describing the 
state of an individual physical system, we do not assign to the Wigner function 
any statistical meaning, and therefore it may safely have negative values: in our 
interpretation, the "negativity problem" is no problem at all. 

Now we want to study the difference between the time evolution of a quantum 
system in the density matrix representation and the time evolution of the corre- 
sponding classical system in phase space. We will use the Wigner representation, 
but the results may be easily carried to position space and momentum space, since 
we will express them in operator form. The quantum evolution, obtained from 
(0), is simply given by: 



dFq 



1 P5 Pd + 1/ (q5 + ^Qd) - ^ (Q5 - ^Qd 



Fq (19) 



To obtain the corresponding classical evolution we start from the continuity equa- 
tion for the classical joint probability density in phase space: 

^ + f^Ji-f^V- = (20) 

at ox m op 

where V is the space derivative of the potential V. Translating ( pO| ) to operator 
form we then obtain 



ih^ = HcFc 



Fc (21) 



Expressions (|l9| ) and (^) are quite similar; indeed, if we expand V (x) in powers 
of X, we can write: 

Hq-Hc = E (Qs) (22) 

The difference between the quantum evolution and the classical evolution involves 
the space derivatives 1/(2"+^) with n > 1. Therefore the quantum evolution and 
the classical evolution are equivalent V = (free particle), V = kx (particle 
moving in a uniform force field) and V = kx^ (linear harmonic oscillator). This is 
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a rather surprising result since we know that the quantum energy spectrum of the 
hnear harmonic oscillator is discrete, and this is a highly non-classical feature; we 
will examine this point in detail in the next section. 

Now we turn to the problem of finding localized solutions in the density matrix 
representation; the existence of such solutions is fundamental for our approach. 
It is well known that to each density matrix satisfying the Schrodinger equation 
(^) we can associate a quasi-classical trajectory, namely the time evolution of the 
mean values {x (t)) and {p{t)) of position and momentum; this is of course a 
consequence of the Ehrenfest theorem: 

{x) = -{p) {p) = -{V'{x)) (23) 

m 

From the trivial relation 

{V'{x)) = V'iix)) + y("+i)((:r)) {{x - (x))") (24) 

n=2 



it then follows that the equations (23) define a classical trajectory only if the 
derivatives V^"^^ vanish for n > 2 (again the free particle and the linear harmonic 
oscillator), while in the general case the trajectory will be quasi-classical if the 
moments ((x — {x))^) are small, i.e. if the wave-packet is well localized in position 
space around its mean value {x (t)). 

In standard quantum mechanics, due to the Heisenberg uncertainty principle, 
wave-packets normally spread out with time; only in some special cases this does 
not happen, the most famous example being the coherent states of the linear 
harmonic oscillator first introduced by Schrodinger in 1926 However, in the 
density matrix representation the Heisenberg uncertainty principle is not true, 
and therefore it should be always possible to find solutions which do not spread 
out with time and remain well localized around their center of mass {x (t)), in the 
sense that their "matter density" (p {xs,X£),t) is significantly different from 

xu=0 

zero only in a small region around xs = {x{t)). In our approach, the existence 
of such solutions is of the utmost importance and allows us to state the following 
"localization postulate" : the only physical solutions of the Schrodinger equation 
(Q) are those which remain well localized around their center of mass {x (t)) in the 
limit t — > iboo; by "physical" we mean "representing individual physical systems", 
as opposed to statistical mixtures or ensembles. In the rest of this paper, we will 
use the definition "quantum trajectories" to label such solutions: thus a quantum 
trajectory is the time evolution of a quantum matrix, as defined at the beginning 
of the present section. 

Clearly, our localization postulate is expressed as a qualitative statement, since 
the concept of "well localized" wave packet is not sharply defined. However, even 
in its rudimentary form, this postulate implies an important difference between 
our approach and the standard interpretation of quantum mechanics: in standard 
QM, immediately before the measure of an observable O (for instance the posi- 
tion), the state of the system is usually a linear superposition of eigenstates of 
the associated hermitian operator; then, when the measure is performed and the 
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result On is obtained, the state suddenly collapses to the eigenvector \0n). This 
wave function collapse has always been an obscure feature of the standard QM 
interpretation: many alternative explanations have been proposed (many worlds 
splitting environment induced decoherence [|^, U], spontaneous localization 
[0, |l^, ...) but until now none of these alternative explanations seems to be 
universally accepted by the academic community. 

On the contrary, in our approach the physical wave-packets are always well 
localized around their center of mass {x (t)) and therefore no wave function collapse 
is needed. When a measure of position is performed, and the result is a certain 
position xq, this means that the wave-packet immediately before the measure was 
already well localized around the point xq: our localization postulate explicitly 
prevents the spreading of the physical wave function cp {x,y,t) over a wide space 
region as t ^ cxd. 

Besides eliminating the need of the wave function collapse, our localization 
postulate implicitly negates the validity of the superposition principle: in our 
approach, a linear combination of physical solutions is not a physical solution 
any more, even if it satisfies the linear motion equation (|^. This is a trivial 
consequence of the fact that in general an arbitrary superposition of localized 
wave-packets is not a localized wave-packet. 

However, a linear superposition of quantum trajectories may have a statistical 
interpretation. Let's consider the set {fx} of all quantum trajectories, where A is 
some appropriate index (for instance, A may be the center of mass and momentum 
at time t = 0, together with some non-classical internal state). Then the integral 



99 (x, y,t)= J f (A) (x, y, t) dA (25) 
where / (A) is a non-negative probability distribution, represents a statistical en- 



semble of particles; in our approach, the solutions defined by (P5[) play the same 
role as the density matrices in standard quantum mechanics. Besides, our ap- 
proach rejects those solutions of the Schrodinger equation (^) which are neither 
quantum trajectories nor positive superpositions of quantum trajectories: since 
they do not represent neither individual particles nor statistical ensembles of par- 
ticles, we conclude that they are just mathematical objects with no physical mean- 
ing. 

If we now apply the Moyal-Wigner transformation to the density matrix ([25|), 
we obtain the Wigner function 



F{xs,ps.t)= J f{X)Fx{xs,Ps,t) dA 



(26) 



where Fx {xs,Ps^ t) are the individual Wigner functions associated to the quantum 
trajectories cpx {x, y, t). Note that both the individual functions Fx and the average 
function F may safely have negative values, since neither represents a probability 
density. The only probability density in ( p6|) is the function / (A), which is non- 
negative by definition. 

Let's now examine the problem of finding the quantum trajectories for a speci- 
fied potential V (x). We will suppose that V (x) vanishes as x — > ±00: in classical 
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terms, this means that the particle is subject to a force only in a limited space 
region, while outside this region it moves freely at constant speed. The first step 
consists in writing the Schrodinger equation (13) in the form 



a2 



m 



V \xs+ ^XD 



V [xs - ^XD 



(27) 



From (|27| ) it is easily seen that, given two arbitrary functions / {xs, t) and g {xD,t), 
we may formally impose the boundary conditions 

^ = f{xs,t) and 
xd=o ox 



ip{xs,XD,t) 



D 



{xs,XD,t) 



xs=0 



g{xD,t) (28) 



and then we may extend (^) to a complete solution ip {xs,X£), t) by means of the 
Taylor expansion 



Lp{xs,XD,t) 



n! 9xn 

71=1 ^ 



X 



a; 15=0 



D 



(29) 

with respect to 



where the derivatives are obtained by differentiating n times (| 
xd and then imposing x/j = 0. 

At first sight, this result is rather paradoxical: it implies that we may choose 
an arbitrary time evolution for the center of mass of our system and then obtain a 
solution of the Schrodinger equation which satisfies this time evolution! However, 
there are at least two reasons why the above procedure does not always produce 
acceptable results: the first reason follows from the fact that an arbitrary solution 
of the Schrodinger equation (0) must be a linear combination of the form 



^{x,y,t) = dki dk2 C{ki,k2) fki{x) /fc2 (?/) e ^ 



(30) 



where (x) are the energy eigenvectors of the original Schrodinger equation (||) 
with associated eigenvalues Ek\ of course, if the energy spectrum is discrete then 
the integrals in (|30| ) must be replaced by summations. Therefore the function 
/ {xs, t) in ( |28| ) cannot be chosen arbitrarily, since it must be equal to (30), in the 
case X = y = xs, for some choice of the coefficients C{ki,k2)- For instance, we 
know that for the linear harmonic oscillator the solutions must be periodic in time, 
i.e. ip {x, y,t) = ip {x, y,t + Tq) where Tq is the revolution time of the oscillator; if 
we choose a function / {xs,t) which is not periodic, then it will be impossible to 
obtain from it a complete solution cp {xs, XDjt) by means of (^), i.e. the Taylor 
series ( p9|) will diverge. 

The second reason is that the solutions obtained from ( pTj ) and (]2^), even if 
they are well defined from a mathematical point of view, will in general violate 
some fundamental physical principle, such as the Ehrenfest theorem (^), the 
energy conservation principle or the unitarity condition 



d_ 

dt 



p{xs,XD,t) 



XD=0 



dxs = 



(31) 



which may be regarded as an expression of the non-relativistic principle of mass 
conservation. Therefore, our next step consists in stating clearly the assumptions 



8 



under which the above cited principles may be derived from the Schrodinger equa- 
tion; we will then require that our quantum trajectories satisfy these assumptions. 
The unitarity property pi]) may be derived from the following condition: 



dip 



dx 



D 







as 



±00 



A sufficient condition to obtain the Ehrenfest theorem (p3[) is 



xs 



dip 



dx 



D 



and 



a; 13=0 



dx\, 







as 



xs 



ibcxD 



and finally the energy conservation principle follows from 

d'^p 



dx'jj 







as 



xs 



±00 



(32) 



(33) 



(34) 



Thus the assumptions (p2[)-(p^) depend only on the field p and its three first 
derivatives with respect to xjj computed at = 0. Besides, from the expressions 
(p^-(p^) we see that the same is true for the three observable quantities, i.e. 
center of mass, momentum and energy (in this case only the first two derivatives 
are involved). If we now set 



(0) 



Xr,=0 



and 



in) 



d'^p 



dx 



D 



n 



1,2,3 



X]j=0 



(35) 



where p^^^ and p^'^'^ depend only on xs and t, we easily derive from the Schrodinger 
equation (|2^) the following relations 





m dp^^^ 


dxs 


^ h dt 




,m dp^^^ 


dxs 


^ h dt 




m dp^'^^ 


dxs 


* n dt 



+ -^V {xs)p 



(0) 



+ 2-^V' ixs)p 
n 



(1) 



(36) 
(37) 
(38) 



From the above considerations we deduce that we do not need to know the com- 
plete quantum trajectory p {xs, XD,t)- all properties of physical interest follow 
from the knowledge of the four functions 99^*) {xs,t), with i = 0, ... ,3. We only 
need to impose that these four functions satisfy the relations (^6|)-(^) together 
with the boundary conditions (^2|)-(^), and that p^^^ (xs^t) be a linear combi- 
nation as in (|30|), with x = y = xs, for some choice of the coefficients C{ki,k2). 
Surprisingly enough, the relations (|36|)-(p8|), obtained from the quantum evolution 
(19), are exactly the same that we would obtain from the classical evolution ([2l]). 
Therefore, if a classical trajectory p^^^ {xs,t) = 5 {xs — Xd (t)) may be expressed 
as a linear combination (|30D, then it can be exactly reproduced in the quantum 
domain: this is more likely to happen in the case of open trajectories, for which 
the quantum energy spectrum is continuous, rather than for closed trajectories. 
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where the quantum energy spectrum is discrete and most classical energy levels 
are forbidden; let's examine these two cases in more detail. 

An open trajectory represents a particle which in the limit t — > ±00 behaves 
like a free particle; at some finite time, the particle interacts with the potential 
V {x) and undergoes a scattering process. In one space dimension, there are only 
two possible outgoing directions, i.e. the particle may be transmitted or reflected; 
this is true both in the classical and in the quantum description. However, in the 
quantum domain we must take into account a new, highly non-classical feature, 
namely the tunnel effect: particles with energy less than the potential peak Vmax 
may cross the barrier, and viceversa particles with energy greater than the bar- 
rier peak may be reflected. Therefore, in the quantum case we have to remove 
the classical restriction that all particles with kinetic energy vnv^ jl > Vmax are 
transmitted while all particles with mv"^ /2 < Vmax are reflected. Thus, knowing 
the energy of a quantum particle is not enough to determine whether it will be 
transmitted or reflected; since we believe that the physical laws are deterministic, 
we deduce that the property of an incoming wave-packet to be transmitted or 
reflected depends also on the form of the wave-packet, not only on its energy: in 
other words, besides the center of mass and momentum, the state of a quantum 
particle must include also some internal, or hidden, degree of freedom even in the 
case of a spinless uncharged particle. 

Furthermore, the existence of non-classical tunnelling particles implies that 
their wave-packet must lose its sharp localization while crossing the barrier region; 
if this were not true, then (p3|) and would yield a quasi-classical trajectory for 
the center of mass {x {t)), thus preventing the appearance of non-classical effects. 
This is the reason why in our localization postulate we explicitly introduced the 
limit t ±00: in general, an open quantum trajectory will be well localized 
when t — > —00, i.e. when it is still far away from the potential region; then, while 
interacting with the external potential, it will spread out to some extent; finally, as 
t +00 and the particle leaves the potential region, it will recover its localization. 
In the next section, we will see a practical application of our approach to the tunnel 
effect. 

Let's now consider briefly the case of closed trajectories. It is natural to assume 
that quantum closed trajectories are periodic, i.e. 93^ {x,y,t + Tci) = ^p^i {x,y,t) 
where Td is the revolution time. The following condition must then be fulfilled 

1 27r 

- (Eki - Ek2) = n— (39) 

for some integer number n, where again Eki, Ek2 are energy eigenvalues of the 
Schrodinger equation for pure states. Since normally the revolution time Td de- 



pends on the energy of the particle, (39) may then be interpreted as a quantization 
condition for the energy spectrum in the quantum domain; an important exception 
is provided by the linear harmonic oscillator, where notoriously the revolution time 
is independent from the energy. In our approach the condition (^) is responsible 
for the existence of discrete energy spectra. 
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3 Applications 
3.1 Free particle 

In the free particle case, the quantum trajectories are simply given by 

(x, y, t) = 5(xs-xo- vo t)ei (40) 

that is, the quantum description is equivalent to the classical description. We 
will briefly show, by means of the Moyal-Wigner transformation, the classical 
probability densities associated to some well known quantum states. 
In the case of a plane wave we have 

V;(x,t) = e* 1^(^-1^*) =^ F{x,p,t) = S{p-po) (41) 

Therefore, the plane wave describes an ensemble of classical particles moving at 
speed vq = ^Po uniformly distributed over the x axis (the normalization is one 
particle per unit of length) . 

In the case of a particle perfectly localized at x = 0, whose initial quantum 
state is given by ■0 (x, 0) = 27rh 5 (x), we have: 

^{x,t) = je'n{^-^*)dp =^ F{x,p,t)=5{x-^t^ (42) 

Thus we obtain an ensemble of classical particles following the trajectories x = vt, 
with V = ^p, uniformly distributed in momentum space (the normalization is one 
particle per unit of momentum). 

Finally, for a gaussian state we have: 



2 p2 (x~p/mtf 



^ix,t)=J^^ e 4Ax^a(t) ^^(^.^p^i) 1 g 2Ax,?^ ' 2Axc? ' (43) 

y ^/2^^Axo a{t) vrft, 

where a{t) = 1+ ^J^^ , t. The marginal deviations for the classical joint dis- 
tribution F{x,p,t) at t = are Axq and Apo = ^^^^ , yielding the minimum 
deviation product allowed by the Heisenberg uncertainty relations. For t>0 Ap 

remains constant, while Ax grows according to the law Ax = y Axq^ + ^^2^^^i 
in agreement with the quantum predictions. 

3.2 The linear harmonic oscillator 

Let's consider the case of the linear harmonic oscillator, defined by the quadratic 
potential 

V (x) = ^ mwgx^ (44) 

In the standard representation of the Schrodinger equation (|^), it is well known 
that the energy spectrum is discrete, and the energy eigenvalues are 

En=(n+^jhuJo n = 0,1,2,... (45) 
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On the contrary, in Section Q we saw that for a quadratic potential the quan- 
tum evolution in the density matrix representation is equivalent to the classical 
evolution: the physical quantum trajectories are then given by 

{x, y, t) = 6[xs- x,i (t) ) ei (46) 
where pd and Xd belong to classical trajectories: 



Pd (t) = ^/2mEo cos {uJot + (/>o) Xci(t) = — \ sin(tJoi + 0o) (47) 

Wo V 

Therefore, all classical energy levels Eq are also acceptable in the quantum domain, 
yielding a continuous energy spectrum; this seems to be a severe contradiction be- 
tween our approach and ordinary quantum mechanics. We will now show that this 
contradiction is produced by an inherent limitation in the formalism of ordinary 
quantum mechanics. 

In the Schrodinger equation for pure states, the hamiltonian operator plays 
two different roles: on one hand it generates the time evolution of the system, on 
the other hand it is the energy operator. Therefore its eigenvectors are at the same 
time the stationary states of the system and the energy dispersion-free states; its 
eigenvalues are at the same time the natural frequencies of the system and the 
possible outcomes of an energy measurement. If the natural frequencies form a 
discrete set, then the energy spectrum must be discrete too. 

On the contrary, in the density matrix representation we have two different 
operators: the hamiltonian H, as expressed by (|3|) or by (^), generates the time 
evolution of the system, and therefore its eigenvalues are the natural frequencies; 
the operator Es = 2m ^(Qs)' which appears in the definition ( [l^ ) of the 
total energy, is the "physical" energy operator, and its eigenvectors are the energy 
dispersion- free states. Since in the Wigner representation the operator Es is a c- 
number, its eigenvalues form a continuous set and include all the classical energy 
levels; this is of course a consequence of the commutation relation [Qs,Ps] = 0, 
and is true for all quantum systems, not only for the linear harmonic oscillator. 
Does this mean that all classical energy levels are always physically realizable 
in the quantum domain? No, because in general the energy eigenvectors do not 
belong to physical quantum trajectories, as defined by our localization postulate 
and by our quantization condition (^9[). However, in the specific case of the linear 
harmonic oscillator, the answer is yes: all classical trajectories are also quantum 
trajectories, and all classical energy levels are physically realizable in the quantum 
domain; indeed, since the revolution time does not depend on the energy, the 
condition ( p9| ) does not impose any restriction on the admissible energy values. 

As for the natural frequencies, they form a discrete spectrum already in the 
classical statistical description of the linear harmonic oscillator: indeed, the classi- 
cal probability density F {x,p, t) is periodic, i.e. F {x,p, t + Tq) = F {x,p, t) where 
To = ^ is the revolution time of the oscillator. Therefore its Fourier expansion 
contains only the integer harmonics of the fundamental frequency ujq; these are ex- 
actly the same observable natural frequencies provided by the standard quantum 
description, i.e. the differences between two energy levels as defined by (p5|). 
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Thus we conclude that standard quantum mechanics predicts the correct nat- 
ural frequencies but does not in general predict the correct energy levels: the 
equivalence between natural frequencies and energy levels is not a physical law, 
rather it is a consequence of an inherent limitation of the standard quantum 
mechanical formalism. If we switch to the density matrix representation, and 
eliminate the requirement that a quantum matrix be a product of pure states 
ip{x,y) = ip {x) ip* (y), then we should obtain the correct quantum energy levels 
by imposing our localization postulate and our quantization condition (|39[). 

The validity of our approach may be confirmed by examining the ground state 
of the linear harmonic oscillator in the usual quantum mechanical description. It 
is well known that this state, corresponding to the energy level E = ^huiQ, is 
described by the gaussian function 

^(x,t)=(^)%-^^'^-^^* (48) 

If we apply to this state the Moyal-Wigner transformation, we obtain again a 
gaussian function 

F{x,p,t) = -^e ^ ° J (49) 

which is always positive and therefore may be interpreted as a joint probability 
distribution in classical phase space: it reproduces the same statistical predictions 



of the quantum state (48) for measurements of position and momentum. The 
Wigner function (49) represents indeed a stationary state, since the right hand 
side does not depend on time; however, it is certainly not an energy eigenvector. 
On the contrary, it represents an ensemble of classical particles moving with all 
possible energies; the energy mean value is (E) = ^hooo as expected, but the 
standard deviation is AE = ^hujo, while we would expect it to vanish for an 
energy eigenvector. This result confirms our belief that the stationary states are 
not necessarily energy dispersion-free: the equivalence between stationary states 
and energy eigenvectors is just a shortcoming of the usual quantum mechanical 
formalism, and disappears when we switch to the density matrix representation. 

As for the standard energy eigenstates with E > ^hujQ, it is well known that 
the corresponding Wigner functions are negative in some phase-space regions. 
Therefore in our approach they have no physical meaning, i.e. they do not represent 
neither individual particles nor statistical ensembles. 

3.3 The delta barrier potential 

Let's now consider a system whose potential is given by: 

V{x) = Vo6{x) (50) 
A set of orthonormal eigenvectors may be chosen as follows: 

fkix) = -^sinkx fi{x) = ^ cos {k\x\ - (f>k) (51) 
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where A; > 0, while (j)^ is defined by 

h^k 2 
The energy eigenvalues associated to the above solutions are simply 



tan = -^tq-t <(j)k < - (52) 



= ^ (53) 

as in the free particle case. A linear combination of the two states fl and /| gives 
the well known state 

r e^^^' + Aij 6-^^=^ if x<0 
[ At e*''^ if X > 

where 

A "i^O A A ^^^^ /^,-\ 

^« = ^OTP = ^OTP ^ (55) 

— mVo ikn — mVo 

Since a term e^^^^ represents a plane wave with momentum p = ±hk, the state 
( |54[ ) is usually interpreted as follows: ^0 e*^^ is the incident wave, moving from 
the left towards the barrier; An e~^^^ is the reflected wave, moving away from 
the barrier with negative momentum p = —hk; At e^^^ is the transmitted wave, 
moving over the barrier with the same positive momentum as the incident wave. 
However, classical mechanics predicts that no particle can be transmitted over an 



infinite potential barrier; therefore, the state (54) is an example of a non-classical 
feature of quantum mechanics, the so called tunnel effect. In the rest of this 
section we will examine the tunnel effect in the density matrix representation. 

Our goal will be to find quantum trajectories representing transmitted and 
reflected particles; these quantum trajectories will be obtained as solutions of the 
simplified equations (|3^)-(|38|) and will satisfy the fundamental physical principles 
defined in Section |2|, i.e. unitarity, energy conservation and Ehrenfest theorem. In 
the case of a reflected particle, there is a very simple solution: we already know 
that the equations (36)-(|38[) are equal to their classical counterparts, therefore we 



will choose the classical solution: 

(xs, t) = 5 {xs - vt) e i-t) + 6ixs + vt) e (t) (56) 



where f > is the speed and 9 is the Heaviside step function. If we apply ( pq ) to 
(p^^\ we obtain for the momentum density V {xs,t) the following expression: 

V {xs, t) = -ihif'^'^^ {xs, t) = mv(5 {xs - vt) 9 {-t) -5{xs + vt) 9 {t) ) (57) 

as expected, i.e. the incident wave-packet has momentum P = mv while the 
reflected wave-packet has momentum P = —mv. Now we integrate (^) to obtain 
(/?(^); the energy density £ {xs,t) is then given by: 

£ {xs, t) = -^ V^'^ i^s, t) + V {xs) v'^^) 

= i mt;2 (5 {xs - vt) 9 {-t) +5{xs + vt) 9 {t) ) (58) 
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The derivation of ( |5q ) is not straightforward, due to the products of distributions 
appearing in the intermediate calculations. However, it may be easily obtained 
by replacing ( [5^ ) with the classical solution associated to the triangular potential 
{x) = ^ (^1 — 0(^1 — -i^^ and then taking the limit as e — > 0; of course, as 
e — > the function (x) approaches the delta barrier potential ([50|). 

Let's now turn to the case of a transmitted particle. As a first attempt, we 
choose the classical solution for a free particle: 

= 5(xs-wt) (59) 

even if we know that it can't be the final solution, since it violates the Ehrenfest 
theorem at t = 0. From ( |36| ) we then obtain 

^^'H^s,t)=i^6{xs-vt) (60) 
n 



which is again equal to the free particle case, since equation ( pjq ) does not in- 
volve the potential V (xs)- The next step would be to integrate ( |37| ) and find an 
expression for 93^^^; this would produce the following result: 

95^'^ (xs, t) = ^ (^-mv^6 {xs - vt) + Vo 6 (xs) 6 (vt) + Vo 6' (vt) 9 {xs)) (61) 

which does not satisfy the condition ( ^3|) at t = and therefore violates the Ehren- 
fest theorem; responsible for this violation is the step function 6 (xs) appearing 
in the r.h.s. of (^). Similar problems arise when we try to integrate (|38|) : the 
resulting expression for ip^^^ ixs,t) does not satisfy condition (^) at t = 0, which 
means that energy conservation is violated too. 

Thus we have shown by direct calculations what we already knew from phys- 
ical considerations: the classical free field solution ( ^9|) may not represent a non- 
classical tunnelling particle. However, since we now know the mathematical source 
of the above violations, we may find a way to circumvent them: we simply must 
add countertems to (|59|), so that the unacceptable terms appearing in the expres- 
sions for (p^"^^ ixs-,t) and (p^^^ {xs,t) are canceled by the new ones. The resulting 
expression is: 



^(°) {xs, t) = 5ixs- vt) - '-^^ 6 (vt) + %^ 5' (vt) (62) 

or (0) gr (0) 

where gi {xs) and (72 (xs) are finite approximations of the Dirac 6 distribution, i.e. 

even functions localized around X5 = 0; and are their n-th derivatives. 

From the fact that gi and g2 are even we deduce that the derivatives g^^ (0) and 

92^^ (0) vanish when n is odd. The expressions for ip^^^ {xs,t) and (p'^'^^ (xs^t) are 
then: 

t) = z7 (5 {xs - vt) + 4!/^ {vt) - 6" {vt)) (63) 

(^., t) = -^(5 {.s -t) - ^ 5" {vt) + ^ .(3) 
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The step function has disappeared from (|6^), so that (X5,t)and ^(2) ^xs,t) 



now satisfy both conditions (|32|) and (33), required respectively by the unitarity 



principle and by the Ehrenfest theorem. As for ip^^^ (xs,t), we do not write its 
final expression, since it does not correspond to an observable physical quantity; 
it is enough to say that it now satisfies the condition (^) required by energy 
conservation. From (|6^ ) and (|6^ we may then obtain the momentum density 
V {xs,t) = —ihif^^^ ixs,t) and the energy density £ ixs,t) = — ^ V'^^^ {xs:t) + 
{xs,t). 

We further note that the wave-packet defined by (^) is sharply localized at 
xs = vt when t ^ 0, while it has a non-vanishing dispersion around xs = 
when t = 0, i.e. when it crosses the barrier. This temporary loss of localization 
was expected, since it has been already recognized in the previous section as a 
necessary condition for satisfying the Ehrenfest theorem in the case of a non- 
classical tunnelling particle. 

From what we have seen until now, we may say that ( |56|) and (|62|) represent 
quantum trajectories associated respectively to reflected and transmitted particles, 
since they satisfy our localization postulate together with the three fundamental 
physical principles. It remains to show that they may be extended to complete 
solutions of the Schrodinger equation, i.e. that they may be expressed as linear 
combinations 

y,W(x5,i)= E /d^i /dA:2C,,(A;i,fc2)/L(x5)/4(x5)e-H(^'=i-^-)* (65) 

for some choice of the coefficients Cij{ki,k2)- We will not demonstrate it 
here; however, in Appendix A we will explicitly build numerical solutions of the 
Schrodinger equation having a very strong resemblance with (^) and (|62|). 

Thus we have shown that in a typical quantum system, exhibiting the highly 
non-classical behaviour known as tunnel effect, it is still possible to find solutions 
which remain sharply localized around their center of mass as t — > ±00; these 
solutions represent wave-packets which are totally reflected or totally transmitted. 
It is our belief that they represent the individual physical particles; the usual 
quantum wave-packets, which split into a reflected part and a transmitted part, 
are just a statistical mixture of our quantum trajectories. Besides, those quantum 
states which may not be written as positive superpositions of quantum trajectories 
must be rejected as physically meaningless. 

As an example, we will examine again the state defined by (|5^. If we consider 
the density matrix (p{x, y) = ijj (x) ip* (y), we may easily compute the expressions 
for if^^^ {xs,t), V {xs,t) and £ {xs,t); in the case X5 > we obtain: 

(^(0)(x5,t) = \At\' (66) 
V{xs,t) =\ATfhk (67) 

£{xs,t) =\At?^ (68) 

The physical interpratation of (|6^)-(|6^) is straightforward: they describe an en- 
semble of transmitted particles uniformly distributed to the right of the barrier 
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(the particle density being \At\'^ per unit of length), moving with speed v = hk/m, 
momentum p = hk and energy E = fi'^k'^ /2m. To the left of the barrier, i.e. in the 
case Xs < 0, we have: 



ip^-' {xs, t) = \Ao\^ + \Ar\'' + 2Re {AoAy^""^^ } (69) 



,(0) 

|2 



V{xs,t) =\Ao\'hk -\AR\'hk (70) 

Ei.s.t) =Mol=^ + I^Hl^^ (Tl) 

The terms involving \Aq\'^ and {ArI"^ represent respectively an ensemble of in- 
coming and of reflected particles, with the right values for momentum and en- 
ergy, l^ol^ and \Ar\'^ being the associated particle densities; besides, the relation 
|At|^ + |^_rP = l^oP ensures that the total number of particles is conserved at the 
two sides of the barrier. However, the last term at the r.h.s. of ( |69|) has no phys- 
ical explanation: it does not contribute to the momentum and energy densities, 
so it seems to describe an ensemble of particles at rest; yet, it has both positive 
and negative values depending on xs, and therefore it may not be interpreted as a 
particle density. For this reason, we are forced to reject the state (54) as physically 



meaningless, since it is not a positive superposition of quantum trajectories. 

On the contrary, starting from our quantum trajectories ( |5^ ) and (^) it is 
very easy to build a statistical ensemble having the same physical properties as 
(|5^, but with no strange unphysical terms. We simply have to write: 

(/p(o) {xs, t) = (\Ar\' v^g) {xs, t-to) + \At\^ {xs, t - to)) V dto (72) 

where of course is the "reflected" quantum trajectory (l56| ) while ipj^^ is the 
"transmitted" quantum trajectory (^). This is another example of the limitations 
inherent to the standard formalism of quantum mechanics, which may be overcome 
by switching to the density matrix representation. 



Finally, we point out that the quantum trajectories (^) and (62) must not be 
considered as the only possible quantum trajectories for the delta barrier potential: 
specifically, it is well known that standard quantum wave-packets impinging onto 
a potential barrier may experiment a shift in time, either a time delay or a time 
advance. For instance Nakazato shows, precisely in the case of the delta 
barrier potential, that a gaussian wave-packet under certain conditions splits into 
a reflected and a transmitted wave-packet which, in the limit t — > oo, are again 
approximately gaussian; both wave-packets are delayed in time with respect to 
the ideal case of an instant transmission or reflection at the barrier. Therefore, 
it is quite possible that quantum trajectories may be found which are shifted in 
time during the interaction with the barrier; this is again a highly non-classical 
feature of tunnelling particles, since of course in the classical case the interaction 
time with a delta barrier vanishes and no time shift may occur. 
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4 Discussion 



In the early days of quantum mechanics, Schrodinger tried to interpret the solu- 
tions of his equation as "matter waves" , i.e. he thought that the individual particles 
could be described by sharply localized wave-packets. Unfortunately, this realistic 
interpretation of the wave- function was soon abandoned, when it became clear that 
wave-packets lose their localization with time and become indefinitely extended as 
t — go; this is clearly incompatible with the fact that particles are always detected 
in small space regions. Then Born proposed his interpretation of the wave-function 
as "probability waves" , and his statistical postulate, incorporated in the standard 
Copenhagen interpretation, has survived until today. The theory built upon this 
postulate, standard quantum mechanics, gives predictions which are in spectacu- 
lar agreement with all experiments ever performed and was successfully extended 
to include relativistic phenomena, such as particle creation/annihilation, and to 
describe all the relevant interactions of the microphysical world. 

Nevertheless, the history of quantum mechanics is not only a history of experi- 
mental successes; it is also a history of never ending debates about its foundations. 
Countless interpretations have been proposed, alternative to the standard Copen- 
hagen interpretation, but none of them has been capable of gaining universal 
acceptance inside the scientific community: the supporters of the various inter- 
pretations still debate throughout the specialized literature, emphasizing both the 
merits of their favourite theories and the flaws of the opponent interpretations. It 
is not my intention here to give a detailed description of the fundamental problems 
which still prevent a completely satisfactory interpretation of quantum mechanics; 
I just want to focus the reader's attention about a few concepts which represent 
the basis on which the present paper is built. 

Realism . A physical system exists as an objective reality, independently from 
the presence of an observer performing measurements on it; therefore a sound 
physical theory should be able to describe the objective properties of the system. 
On the contrary, standard quantum mechanics only provides statistical predic- 
tions about the possible outcomes of experiments, denying the existence of ob- 
jective properties prior to measurement: the famous Schrodinger cat is described 
by a superposition of "dead" state and "alive" state, until somebody opens the 
box and the cat suddenly collapses to either "dead" or "alive". Transposing the 
paradox to the description of a single particle, we have a wave-function which may 
be extended over an arbitrarily large space region; then the particle's position is 
measured and the wave-function suddenly collapses to the state |xo), where xq is 
the measurement's result. A much more realistic view would be to think that im- 
mediately before the measurement the particle is already localized at the position 
xq, this position being an objective property of the particle: the measurement 
simply brings to light this property already possessed by the particle. I am aware 
that such way of thinking would be labeled as "naif realism" by the majority of 
physicists; nevertheless, I am still convinced that it is the only acceptable starting 
point for building a sound physical theory. 

Determinism . Probabilities are a very useful tool for handling systems whose 
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state, for whatever reason, is not completely known: they allow us to make calcu- 
lations and to obtain interesting results even if we are partially ignorant about the 
system's state. Standard quantum mechanics has radically changed this original 
meaning of the probability concept: on one hand the wave- function is assumed 
to represent a complete description of an individual physical system, on the other 
hand the Born postulate implies that we cannot predict with absolute certainty 
the result of a measurement performed on a system, even if its wave-function is 
perfectly known. Thus the concept of uncertainty becomes a fundamental princi- 
ple of the theory, and the Heisenberg uncertainty relations play a central role in 
the so-called "quantum revolution" . 

For many years it was generally believed that a fully deterministic theory re- 
producing the same statistical predictions of quantum mechanics could not exist; 
Von Neumann's proof of the impossibility of "hidden variables" theories was con- 
sidered as the last word about the subject. However, in 1952 Bohm Q was able to 
build exactly this kind of theory: in Bohmian mechanics, the state of a quantum 
particle is given by its wave-function ip (x) and by its position Q. The time evolu- 
tion of the wave-function follows the usual Schrodinger equation, while the time 
evolution of the particle's position is obtained from a new equation, satisfying the 
following remarkable property: if we consider at time t = a statistical ensemble 
of particles, all described by the same wave-function ip^XjO), whose positions Q 
are distributed with a probability density p{x,0) = \ip{x, 0)1"^, then the equality 
p{x,t) = \ip {x,t)\^ will remain true for all times t > 0, thus reproducing the same 
probabilities of standard quantum mechanics for measurements of position. Hence 
Bohmian mechanics is at the same time realistic, the position Q being considered 
as an objective property of the particle, and deterministic, since the result of a 
measurement performed on the system is uniquely determined by the system's 
state (here the assumption is made that only those observables which may be 
expressed in terms of position measurements have physical meaning). 

Of course, the existence of at least one "hidden variables" theory equivalent 
to standard quantum mechanics is enough to disprove all claims about the impos- 
sibility of such theories. In the present paper, as in the previous one, I strongly 
support the view that the fundamental laws of nature must be deterministic: if 
there is uncertainty about the results of measurements performed on a physical 
system, this simply means that the state of the system is not entirely known. For 
instance, let's consider two particles impinging onto the same potential barrier, 
both starting with the same initial position and momentum; if the first particle 
is reflected while the second one is transmitted, this is not a consequence of some 
fundamental uncertainty in the laws of nature: on the contrary, this means that 
position and momentum are not a complete description of the particle's state. 

Nevertheless, there is one serious problem that must be faced by every de- 
terministic theory trying to reproduce the statistical predictions of quantum me- 
chanics: of all possible probability distributions over the system's state, why only 
a few may be observed in real experiments, while the majority of them seems to 
be forbidden by nature? In the case of tunnelling, we can imagine a statistical 
ensemble of particles all being reflected (or transmitted); however, the experimen- 
tal evidence tells us that the real ensembles are always formed by both reflected 
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and transmitted particles, and their relative weight in the mixture depends on the 
energy of the ensemble in a well defined manner. In Bohmian mechanics, we could 
think of an ensemble of particles having all the same initial wave- function {x,0) 
and the same initial position Qq; this ensemble would have a dispersion- free posi- 
tion at all future times, thus contradicting the quantum principle that a particle's 
trajectory may not be perfectly known at all times. 

To solve this problem, the supporters of "hidden variables" theories usually 
propose to complete the deterministic equations of motion with some kind of 
statistical postulate, in analogy with the classical case where the thermodynamic 
equilibrium hypothesis (Gibbs postulate) allows to obtain the laws of statistical 
mechanics from the deterministic laws of newtonian mechanics. In the case of 
Bohmian mechanics, this statistical postulate is the already mentioned hypothesis 
that p{x,0) = \ip {x, 0)1^ at some initial time t = 0; in the literature, this postulate 
is usually referred to as the "quantum equilibrium hypothesis" 

Simplicity. The history of physics supports the belief that the fundamental 
laws of nature should be as simple as possible; even if this is not a clearly defined 
principle, it generally means that the fundamental objects of a theory, together 
with its fundamental equations, should be reduced to the smallest possible num- 
ber. For instance, the theory of special relativity unified the electric field and the 
magnetic field into a single object, the electromagnetic field, at the same time 
reducing the four Maxwell equations to a single covariant expression; besides, the 
elimination of absolute space and time brought as a consequence that the laws 
of physics must be the same in all inertial reference frames (Lorentz invariance). 
Then came the general theory of relativity, which further simplified the defini- 
tion of reference frame by eliminating the priviliged role of inertial frames and by 
establishing the equivalence between accelerated frames and gravitational fields; 
as a consequence, the laws of physics must be the same in all reference frames, 
inertial or not (general covariance). Besides, general relativity provides an elegant 
explanation of the gravitational forces as effect of space-time curvature and clar- 
ifies the status of space-time, which becomes a true physical object, not only an 
arena where physical events take place; it is generally believed that the Einstein's 
equations for the gravitational field are one of the highest peaks reached by the 
physical science in its search for simplicity. 

Unfortunately, the same level of simplicity has never been reached by quantum 
mechanics and quantum field theory; on the contrary, there are many controversial 
questions, ranging from the interpretational problems to the divergencies which 
plague quantum field theory and which still prevent the unification with gravity. 
Here I just want to concentrate about one point, which is relevant to the approach 
that I am proposing in the present paper: the so called wave-particle duality. It 
is well known that particles manifest wave-like behaviours in some experimental 
situations: for instance the interference fringes obtained in the famous two-slit 
experiment appear to be a consequence of the linear superposition of two waves, 
adding up where the phase is the same and canceling each other where they have 
opposite phases. At the same time, the individual particles are always detected 
as (almost) point-like objects, i.e. each of them leaves a well localized spot on the 
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screen: the interference fringes may be seen only after a big number of particles has 
been detected. This wave-particle duality is not explained by standard quantum 
mechanics, rather it is considered as a postulate, not requiring further explana- 
tions; embarrassing questions such as which slit did the particle pass through are 
rejected as meaningless, since in standard quantum mechanics there is no room 
for such an old-fashioned concept as the trajectory of a particle. 

On the contrary, Bohmian mechanics provides a clear explanation of wave- 
particle duality: in the case of the two-slit experiment the wave function passes 
through both slits, while the particle trajectory, expressed by the function Q (t), 
passes through but one of the slits. However, the point that I want to underline 
here is that this explanation does not meet the simplicity criterium defined above: 
instead of unifying the two concepts (particles and waves) into one single object, 
capable of manifesting both behaviours, the duality is solved by stating that par- 
ticles and waves exist as two separate entities: the number of fundamental objects 
of the theory, together with its fundamental equations, is increased with respect 
to standard quantum mechanics, while the predictive power of the theory remains 
the same. This is probably the main reason why Bohmian mechanics has been 
considered by the majority of physicists more like a mathematical curiosity, rather 
than a serious alternative to standard quantum mechanics. 

As a counteraxample, I will now briefly outline how the wave-particle duality 
is solved by a recently proposed theory, namely Hasselmann's "metron model" 
In this theory, the only fundamental object is the metric tensor g^M defined over a 
higher- (eight- or nine- )dimensional space, where the first four dimensions define 
the ordinary space-time, while the last four or five are extra-space dimensions. 
The only fundamental equation is the Einstein's gravitational field equation 

Rlm = (73) 

where Rlm is the Ricci curvature tensor. The author shows that the simple equa- 
tion ( |73|) has "an extremely rich nonlinear structure which encompasses all the 
principal interactions of quantum field theory and can be used as the founda- 
tion of a unified deterministic theory of fields and particles". Specifically, it is 
postulated that soliton-type solutions exist, named "metrons", consisting of "a 
localized, strongly non linear core and a set of linear far fields ... The core is the 
origin of the corpuscular properties of matter, while the ... far fields give rise to 
the wave-like interference phenomena". Thus the metron model unifies gravity 
with the other forces of nature and solves the wave-particle duality, and all this 
is achieved starting from equation ([73[), which is even simpler than the original 
Einstein's field equations: indeed the energy-momentum tensor, which appears as 
an external source term in the original field equations for gravity, is not present in 
the fundamental equation of the metron model; on the contrary, the author shows 
that the standard energy-momentum tensor arises from the contraction of the 
extra-space components of the Riemann curvature tensor. Hasselmann's model 
is a striking example of a theory which meets our simplicity criterium; besides, 
it has more than one point in common with the approach that I propose in the 
present paper. 
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So far we have identified three important guidehnes for building a sound phys- 
ical theory: realism, determinism, simplicity; now let's see how the approach 
presented in this paper tries to meet these criteria. Basically, what I propose is to 
revive the original Schrodinger's idea of the wave-function as "matter waves" . The 
main obstacle to this idea, i.e. the fact that wave-packets spread out with time, is 
overcome by choosing as the fundamental equation of the theory the Schrodinger 
equation for density matrices, instead of the one for pure states on which standard 
quantum mechanics is based. In Section |2| we have seen that in the density matrix 
representation the Heisenberg uncertainty principle is not true, i.e. there are fields 
perfectly localized both in position space and in momentum space; as a conse- 
quence, wave-packets do not necessarily spread out with time. We have further 
postulated that the only physical solutions (here "physical" means "representing 
individual particles") are those which remain well localized as t — > ±00, labeled 
"quantum trajectories" . In Section ^ we have seen that quantum trajectories exist 
in three important cases: the free particle, the linear harmonic oscillator and the 
delta barrier potential; I am firmly convinced that they may be proved to exist in 
all cases of physical interest. 

Our approach is deterministic: if we know the state if (x, y, 0) of the particle at 
i = 0, then the Schrodinger equation (^) determines the state ip (x, y, t) at all future 
times; besides, the values of the observable quantities (center of mass, momentum 
and energy) are uniquely determined from the state by means of equations (13)- 
(15). Note that the particle's position is not a well defined concept for a wave- 
packet, and therefore we replace it with its center of mass; however, since we 
postulate that the real wave-packets remain sharply localized as t 00, the 
wave-packet's center of mass is practically indistinguishable from the position of 
a pointlike particle. 

Our approach is also realistic: no "Schrodinger cat" state is allowed and no 
wave-function collapse is needed; the position of the particle, i.e. the center of mass 
of a sharply localized wave-packet, is known prior to measurement and therefore 
may be considered as an objective property of the particle. Of course, since the 
state (/9 (x, y) depends upon two position coordinates, our approach postulates the 
existence of an objective extra-space coordinate: in addition to the "physical" 
position xs = {x + y)/2, we have the "auxiliary" position coordinate xd = x — y, 
having observable effects only around the point xd = 0. The idea of extra space- 
time dimensions has a long lasting tradition in physics, which was initiated by 
the original Kaluza-Klein theories ||l^, ^ and has been recently revived by the 
already cited Hasselmann's metron model. 

As for simplicity, we have only one fundamental object, the field ip {x,y), and 
one fundamental equation, the Schrodinger equation (^. The wave-particle du- 
ality is solved by postulating that the physical wave-packets are sharply localized 
as t — > ±00, i.e. when the particles are produced and detected: this explains why 
the individual particles are always experienced as (almost) pointlike objects. At 
the same time, the wave-packets may spread out to some extent while interact- 
ing with the potential V {x); this temporary loss of localization is responsible for 
non-classical phenomena such as tunnelling or interference, i.e. for the wave-like 
behaviour of quantum particles. We may easily imagine the form of a quantum tra- 
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jectory in the case of the two-sHt experiment: it starts as a locahzed wave-packet 
as t — > — oo, then when approaching the potential region it spUts into two separate 
wave-packets; each wave-packet passes through one of the shts; then, while leaving 
the potential region, the two wave-packets merge again and as t — > +00 a single 
localized wave-packet is detected, producing a pointlike spot on the screen. There 
is an important difference between our approach and standard QM concerning 
the concept of interference: in standard QM, interference takes place between two 
plane (or spherical) waves defined over wide space regions; in our approach, inter- 
ference is localized in space and time, i.e. it is confined to the potential region and 
to the time interval when the physical wave-packet interacts with the potential. 
Outside this space-time region our wave-packets move along classical trajectories, 
behaving as classical free particles; however, the interaction with the potential 
bends the outgoing trajectories with respect to the classical case, producing sta- 
tistical effects which may not be explained classically, like the interference fringes 
in the two-slit experiment. 

The superposition principle of standard QM does not apply to our quantum 
trajectories, i.e. a linear combination of quantum trajectories does not in general 
represent an individual particle. However, it may be accepted in its statistical 
meaning: a positive superposition of quantum trajectories may be taken to repre- 
sent a statistical ensemble (or mixture) of particles. This statistical generalization, 
of course, is not enough to reproduce the predictions of standard QM; we still need 
some kind of statistical postulate (similar to the quantum equilibrium hypotesis 
in Bohmian mechanics) to explain, for instance, the relative frequency of trans- 
mitted and reflected particles in the case of tunnelling, or the exact form of the 
interference fringes in the case of the two-slit experiment. This remains an open 
issue in our approach; nevertheless, it seems to me that the solution to this prob- 
lem should be made easier by the fact that our fundamental equation, both for 
individual particles and for statistical ensembles, is exactly the same equation 
which describes the time evolution of density matrices in standard QM: we did 
not change the equation, we only adopted a different criterium to decide which 
solutions are physically meaningful and which are not. 

There is a second open issue in our approach: the definition of quantum tra- 
jectory is slightly vague, since we did not specify exactly what we mean by "well 
localized wave-packet" ; should it be a Dirac delta function, with vanishing disper- 
sion, or is it allowed to have a finite extension (for instance a sharp gaussian)? 
Moreover, the Schrodinger equation (^) admits solutions which do not have physi- 
cal meaning, i.e. they are neither quantum trajectories, nor positive superpositions 
of quantum trajectories; we had to introduce a specific postulate to define which 
are the physical solutions, but clearly a truly fundamental theory should contain in 
itself a mechanism for selecting the physical solutions and discarding the unphys- 
ical ones. The origin of this problem is evident: it is rooted in the linearity of the 
Schrodinger equation. It is linearity which allows, by means of the superposition 
principle, to obtain finite dispersion wave-packets starting from delta functions; 
it is again linearity which allows both positive and negative superpositions, thus 
producing solutions for which no physical interpretation is possible, not even at 
the statistical level. 
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Linearity is a much celebrated feature of standard quantum mechanics: it is 
often claimed that without the superposition principle the wave-like features of 
elementary particles could not be explained. The mathematical apparatus of the 
theory is totally based on linear objects (Hilbert spaces, matrix operators, com- 
mutation relations, . . . ); since this linear apparatus happens to work pretty well in 
the non-relativistic case, it has been extended with few modifications to describe 
relativistic phenomena. Even if the experimental successes of quantum field theory 
may not be denied, I do not share all this enthusiasm about linearity; on the con- 
trary my opinion is that linearity is responsible for some of the problems which 
plague quantum mechanics and quantum field theory. At the interpretational 
level, it is the main origin of confusion about the meaning of the wave-function: 
does it represent an individual particle or a statistical ensemble of particles? Of 
course this question may not find an answer inside a linear theory, where both 
quantum trajectories (representing individual particles) and positive superposi- 
tions of quantum trajectories (representing statistical mixtures) are solutions of 
the fundamental equation of the theory! At the formal level, I am convinced that 
the divergences encountered in quantum field theory, with the consequent need 
of cumbersome renormalization procedures (not always effective, as in the case 
of gravity), are originated by the attempt of forcing a linear structure upon a 
phenomenology which is essentially non-linear. 

As a consequence, I don't believe that an attempt to refine our definition of 
quantum trajectory would be much useful, if confined to the linear Schrodinger 
equation. A much more useful effort would be trying to add non-linear terms to 
the equation, with the goal of eliminating the unphysical solutions: the only stable 
solutions of the new equation should be those which represent individual particles. 
The premises for a non-linear modification of the Schrodinger equation have been 
studied by many authors (see for instance [0, ^ ^, |l^]); these explorations have 
been motivated mainly by the general observation that "all linear equations de- 
scribing the evolution of physical systems are known to be approximations of some 
nonlinear theories, with only one notable exception of the Schrodinger equation" 
[||. This is undoubtedly a well founded observation; unfortunately, our motivation 
is much more specific, and the existing literature seems to be of little help in our 
case, especially because the starting point of the above cited authors is usually 
the Schrodinger equation for pure states, not the one for density matrices. Be- 
sides, we must remember that the non-relativistic Schrodinger equation is just an 
approximation for some relativistic equation, the Dirac equation for fermions or 
the Klein-Gordon equation for bosons; it seems to me that the search for a fun- 
damental equation should start already at the relativistic level, i.e. the non-linear 
modifications should be applied directly to the Dirac or Klein-Gordon equations. 
An example of a non-linear modification of the Dirac equation has already been 
proposed in my first paper [18|, where the non-linearity was introduced by coupling 
the Dirac field with the classical electromagnetic field. 

The attentive reader may have noticed that our main effort has been devoted 
to the subject of open quantum trajectories: closed trajectories have been intro- 
duced only briefiy at the end of Section ^ and examined in detail only for the 
linear harmonic oscillator (Section |3|). The main conclusion of the present paper 
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about closed quantum trajectories is the distinction between physical energy levels 
and natural frequencies, two concepts which in standard quantum mechanics are 
taken to coincide; in our approach, the natural frequencies are the eigenvalues of 
the hamiltonian operator H, defined by (|3|), while the energy levels are the eigen- 
values of the physical energy operator Es = 2^ P5 + ^(Qs)- However, since the 
spectrum of the operator Eg is continuum, thus including all classical energy lev- 
els, the fact that the outcomes of energy measurements are quantized is explained 
by our approach in a different way than in standard QM, namely by means of the 
quantization condition (^). Of course this subject needs further investigation; 
here I just want to point out a simple observation, which may have some inter- 
esting physical consequences. The observation is, trivially, that nobody has ever 
seen what happens inside an isolated atom: what we see is how the atom reacts 
to some external excitation, for instance the freqeuncy of the light emitted by the 
atom after it has been excited by the collision with an accelerated electron; the 
hypotesis that the emitted light frequencies correspond to jumps of the internal 
electrons between two different energy levels, as defined by standard QM, may not 
be experimentally tested. It is the only possible explanation if we confine ourselves 
to the Schrodinger equation for pure states; however, in the density matrix repre- 
sentation there may be other explanations, in terms of closed quantum trajectories 
whose total energy does not necessarily coincide with the energy levels of standard 
QM. In the end, the only observable effect is that the emitted light frequencies 
correspond to the natural frequencies of the system under examination, and these 
are the same both in standard QM and in the density matrix representation. 

The present paper has been entirely devoted to the case of a single particle; 
in Appendix B we will briefly sketch a possible way of dealing with many-particle 
systems in the density matrix representation. Here I just want to point out that 
every realistic interpretation of a many-particle quantum system must face the well 
known dilemma of quantum non-locality: the result of a given measure performed 
on one particle seems to be influenced by what kind of measurement is performed 
on a second particle, distant from the first. However, this weird feature of quan- 
tum mechanics is not fully established: there are controversies still open both at 
the theoretical and at the experimental level. On one hand, there are authors 
[ pC| , ^ who believe that the quantum probabilities should be considered as condi- 
tional probabilities and therefore they have no right to enter the Bell inequalities; 
on the other hand there are still doubts about the experimantal violation of Bell 
inequalities, due to the well known loopholes Therefore the impossibility 

of reproducing the experimental predictions of standard QM (and the observed 
experimental results) by means of local realistic theories has not yet been proved 
beyond doubt, and maybe never will be; my attitude about this point is not dif- 
ferent from the one I have taken in the single particle case, where those quantum 
states which could not be expressed as positive superpositions of quantum tra- 
jectories were rejected as physically meaningless: in the many-particle case, I am 
inclined to reject as physically meaningless those quantum states whose statistics 
may not be reproduced by local realistic theories. 
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A Numerical quantum trajectories for the delta bar- 
rier potential 

As a first step, we normalize tlie Schrodinger equation by imposing m = h = VQ = l; 
this choice does not result in a loss of generality, since we may get back the original 
Schrodinger equation through a suitable rescaling of the coordinates x, t and of 
the field (p. As a second step, we must choose a finite set of eigenvectors out 
of (|5l|), since it is impossible to deal numerically with an infinite continuum of 
eigenvectors; therefore we will consider only a limited space region —L < x < L 
and only those eigenvectors / (x) which behave smoothly at the borders of such 
region, i.e. / (L) = / {—L) and /' (L) = f {—L). Our set of eigenvectors will then 
be defined by: 



fn {x) = cos {kn \x\ — (j)n) for n even 
fn (x) = sin knX for n odd 



For n odd, kn is simply given by 



kr, L = — TT 



while for n even it may be obtained by solving the coupled equations: 



tan ( 



kn L 



1 



0<(J)n< 



TT 



For our computations, we will choose L = 100. 

Our numerical quantum trajectories will then have the form 

N N 

ip{x,y,t) = J2 Cijipi,{x,y,t) = Ci,fiix)fj{y)e'^^^~-^^' 

i,j=0 i,j=0 



(74) 
(75) 

(76) 

(77) 
(78) 



(79) 



where uji = kf/2; our goal will be to find coefficients Cij such that ( |79D provides 
good approximations for the simplified solutions (^) and (62), which represent 
reflected and transmitted particles. To accomplish this task, we first associate to 
every quantum trajectory ip {xs,XD,t) the three functions 



p{xs,t) 
Pixs,t) 
e{xs,t) 



ip{xs,XD,t) 



d 



-i- ip{xs,XD,t) 

dXD 



a; 15=0 



1 52 1 52 \ , 



(80) 
(81) 

(82) 



X_n=0 
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which represent respectively the mass density, the momentum density and the 



energy density; the reader may notice that the definition (82) for the energy density 
differs from the usual one obtained from (|T5|). The reason is that, while the total 
energy is the same, the choice made in ( |82|) eliminates from the beginning the 
possible appearance of delta functions at = 0. 
Then we define the scalar product 

{ipiip2) = f dt f dxs (wq pi{xs,t) p2ixs,t) + 
J — T -J — L 

+ wi pi{xs, t) p2{xs, t) + W2 £i(x5, t) e2{xs,t) ) (83) 

where wq, wi and W2 are appropriate weights, while [— T,T] is a time interval 
to be defined. To obtain our numerical quantum trajectories we will project, by 
means of the scalar product (|83|), the simplified solutions (5£) and (B^) over the 



linear space defined by (|79|). For instance, if (pR is the solution associated to a 
reflected particle as defined by (p6[)-(p^), the coefficients Cij will be obtained by 
solving the system of {N + 1)^ linear equations 

N 

( ^ki ^ij) = {^ki ^r) k,l = 0,...,N (84) 

i,j=o 

which requires the inversion of a {N + 1)-^ x (A^ + 1)-^ square matrix. The choice 
of is then dictated by the computational resources of the machine used for 
performing this matrix inversion; in our case we will set = 100. 

For both reflected and transmitted particles we will choose the speed v = 1; 
in the standard quantum mechanical treatment this is the speed for which the 
transmission probability and the reflection probability are the same. Besides, due 
to numerical considerations, we will not use delta functions in the definitions of 
(pR and (fT- On the contrary, we will use the smooth function 

which in the limit Ax — > behaves like a Dirac delta function. In (|85|), 9 is 
the Heaviside step function, while Ax is a finite dispersion to be defined; for our 
computations we will use Ax = 3. Besides, we are now able to choose the time T 
appearing in the definition ( |83|) of our scalar product: we will set T = L — ^Ax, 
so that the wave-packet / (x — vt) is entirely included in the space region [—L , L] 
when t G [-T ,T]. 

Now we are ready to show the results of our computations. In the case of a 
reflected particle, with a suitable choice of the coefficients wq, wi and W2, we obtain 
the quantum trajectory described by figure 1. In figure l.a, we show the wave- 
packet representing the mass density p (xs) at different times t = ±80, ±40, 0. In 
figures l.b and l.c, we show the wave-packets associated to the momentum and 
energy densities at the same times; the momentum density at t = may not be 
displayed, since it vanishes identically for all xg. Finally, in figures l.d and l.e 
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we show the mean value xm and the standard deviation ax of the wave-packet as 
a function of time; these two quantities have been obtained using as weight the 
function p^, for instance the mean value is given by 

dx 

XM = +L 4. . , (86) 

This choice has been made to reduce the effect of the low amplitude noise appearing 
in figure l.a. From figure 1, it is clear that our quantum trajectory is a very good 
approximation of the corresponding classical trajectory. Specifically, the standard 
deviation is approxmately constant when t is different from zero, but is reduced 
for t = 0, due to the overlapping of the incoming part and the reflected part of 
the wave-packet; this is qualitatively the same behaviour that may be seen in the 
classical case. 

Turning now to the case of a transmitted particle, we make a slight change to 
the definition of our scalar product (|3[) and multiply the integrand term by the 
factor 

p{xs,t) = i-e{^Ax-\xs\) e(^Ax-\t\^ (87) 

that is, we eliminate from the integral the region defined by xs G [— vrAx , vrAx ] 
and t £ [—^Ax , §Ax] . The reason is that we do not want to force any form to 
the wave-packet while it is tunnelling through the barrier; after all, this is a non- 
classical behaviour, and we do not have any a priori knowledge of the deformations 
produced on the wave-packet by this non-classical interaction. With a suitable 
choice of the coefficients wq, wi and W2, we then obtain the quantum trajectory 
described by figure 2. The time evolution of the mean value xm (figure 2.d) is not 
different from the case of a classical free particle; however, from figure 2. a we see 
that at t = the wave-packet is heavily deformed, while from figure 2.e we see 
that the standard deviation increases during the interaction time, as expected. 
Thus we have found two complete solutions of the Schrodinger equation which 



approximate pretty well the simplified solutions (56) and (32). Our numerical 
solutions have some limitations: their range of validity is restricted to the space- 
time region defined by xs G [—L,L] and t G [— T,T], while outside this region 
the wave-packet looses abruptly its localization; besides, the localization of our 
wave-packets is not so sharp even inside the region of validity. However, it should 
be clear that these limitations arise from the fact that we are working with a finite 
set of eigenvectors: if we were able to perform our calculations in the limit N ^ 00, 
the validity of our quantum trajectories could be extended to larger space-time 
regions and to sharper wave-packets, thus approaching the limit L — > 00, T — > 00 
and Ax 0. 
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B Brief introduction to the many-particle case 



In standard QM, a system of N particles is described by a wave-function 
^p (xi, . . . ,xn) and the fact that ip is defined over a configuration space, whose 
dimension depends on the particles' number, is a further obstacle to a realistic in- 
terpretation of the wave-function; therefore this undesired feature must disappear 
in the density matrix representation. To be simple, let's consider a system of two 
particles moving in one space dimension under the effect of an external potential 
V (xi, X2); to avoid possible complications which may arise in the case of identical 
particles, we will further suppose that the two particles have different masses mi 
and 7712- The Schrodinger equation for the density matrix if (xi, X2, yi, 1/2) is then: 

In our approach, the solutions of equation (|8^) are taken to represent individual 
physical systems instead of statistical mixtures. To eliminate the dependence of 
if upon the configuration variables (xi,yi) and (x2,2/2), we would like to impose 
the following separability condition: 

V9(xi,X2,yi,2/2,t) = (xi,yi,t) ip2 (x2,y2,t) (89) 

so that the two functions ipi and 992 would describe separately the time evolution 
of the two particles and could be thought to depend upon just one coordinate 
pair (x,y), as in the single particle case. Unfortunately, the condition ( ^9|) is in 
general too strong, and may be satisfied only in some special cases, for instance 
when V (xi,X2) = Vi (xi) + V2 (X2). However, since we know that the observable 
quantities depend only upon the values of 93 in the region where xi ^ yi and 
X2 ~ y2j we may impose a weaker separability condition: therefore we first define 

Xis = ^{xi + yi) XiD = Xi-yi i = l,2 (90) 



and then require that ( p9| ) be satisfied only far small values of xid and X2_d. 
For xiD = X2D = we then obtain the condition 

^ ^ = Pi{xis,t) p2{x2S,t) (91) 

while at first order in xid and X2d we obtain: 

dtp 



ih 



dxiD 
dip 



dx 



= Viixis,t)p2ix2s,t) (92) 

XiD=X2D=0 

= Pi{xis,t)V2{x2s,t) (93) 

XiD=X2D=0 

If we now extend the definitions (|l^) and ([T^ of center of mass and momentum 
to the two-particle case, we easily obtain 

Q = Qi + Q2 = J xs pi {xs) dxs + J XSP2 (xs) dx5 (94) 
P = Pi + P2 = JVi (xs) dxs + JV2 {xs) dxs (95) 
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where we supposed that both pi and p2 satisfy the unitarity condition 
/ Pi {xs) dx5 = 1. As for the energy, it is clearly impossible to define two sepa- 
rate energy densities, since the potential energy depends on the position of both 
particles; this was true already in the classical case. 

At first order in xid and X2Di the Schrodinger equation (|89|) may be decoupled 
in the two separate equations: 

which are the natural extensions of the condition ( |3^ ) to the two-particle case. 
This decoupling of the Schrodinger equation confirms that solutions satisfying to 
both (^) and (p2|)-(p^) may indeed exist. 

At this point it is clear what should be the definition of quantum trajecto- 
ries in the two-particle case: we will choose those solutions which satisfy both 
( |9TD and (^2|)-(^) and which remain well localized around their center of mass 
as t — > ±00; this localization condition applies separately to both functions pi 
and p2- These quantum trajectories represent individual physical systems, while 
positive superpositions of quantum trajectories represent statistical ensembles; in 
our approach, all other solutions are devoid of physical meaning. Of course, to 
be physically acceptable, our quantum trajectories must also satisfy the three 
fundamental principles, i.e. unitarity, Ehrenfest theorem and energy conservation. 
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